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1.0  iNnooocnoN 

This  report  describes  progress  during  Che  first  year  of  a  research 
program  on  the  Fundamental  Processes  in  Partially  Ionized  Plasmas 
conducted  in  the  High  Temperature  Gasdynamics  Laboratory  at  Stanford 
University.  This  research  is  supported  by  a  grant  from  Che  Air  Force 
Office  of  Scientific  Research  (AFOSR-83-0108)  and  is  conducted  under  the 
direction  of  Professors  Charles  H.  Kruger,  Morton  Micchner,  and  Sidney 
Self.  Three  Ph.D.  candidates  are  currently  conducting  their  doctoral 
research  under  this  program. 

Several  space  power  and  propulsion  systems  of  potential  long-range 
interest  to  Che  Air  Force  involve  partially  ionized  plasmas.  Such 
systems  include  MPD  thrusters,  both  open  and  closed  cycle  MHD  power,  and 
thermionic  energy  conversion.  Although  Che  specific  configurations,  the 
exact  operating  conditions,  and  which  of  Che  competing  systems  will 
prove  Co  be  most  useful  in  the  long  term  remain  to  be  established,  it  is 
important  at  this  time  Co  provide  a  broad  fundamental  research  base  in 
support  of  development  activity.  In  particular,  there  are  a  number  of 
key  issues  regarding  Che  properties  and  discharge  behavior  of  partially 
ionized  plasmas  and  Che  interaction  of  discharges  with  fluid  dynamics 
Chat  need  to  be  understood  before  Che  potential  and  llmltaClons  of 
competing  systems  can  be  fully  evaluated.  In  addition,  it  is  Important 
that  outstanding  young  applied  scientists  be  educated  in  these  areas. 

The  present  research  on  partially  ionized  plasmas  is  discipline 
rather  chan  device  oriented  and  is  currently  focused  on  three  major 
areas : 

1.  Recombination  in  molecular  plasmas 

2.  Discharge  effects:  plasma  electrode  interaction 

3.  Interaction  of  discharges  and  fluid  dynamics 


In  addition,  each  of  these  areas  involves  the  development  and 
application  of  modern  plasma  diagnostic  techniques. 


2.0  PBOJBCT  SOHMARIBS 

Included  in  chls  section  are  suannaries  of  progress  in  each  of  three 
project  areas.  Each  project  suomary  contains  the  following  subsections: 
(a)  Introduction;  (b)  Scientific  Merit;  and  (c)  Status  Report.  Addi¬ 
tional  descriptions  isay  be  found  in  the  publications  listed  in  Section 


2.1  Secoabinatloa  la  Molecular  Plasaas 


Introduction 

The  equilibrium  thermodynamic  properties  and  the  quasl-equlllbrlum 
transport  properties  of  partially  ionized  plasmas  have  been  fairly  well 
worked  out,  at  least  with  regard  to  the  theoretical  aspects,  and  many 
Important  results  are  given  in  the  book  Partially  Ionized  Gases  [1.1]. 
There  remain,  however,  a  number  of  areas  where  plasma  properties  are  not 
now  adequately  understood,  particularly  with  regard  to  nonequilibrium 
plasmas.  Since  a  knowledge  of  plasma  properties  is  fundamental  to  any 
system  description,  it  is  Important  that  this  area  receive  appropriate 
attention. 

The  electron  density  is  often  the  single  most  important  thermo¬ 
dynamic  property;  in  many  nonequilibrium  situations  it  is  determined  by 
finite  ionization  and  recombination  rates.  These  rates  in  turn  depend  on 
Che  constituents  of  Che  gas  and  the  form  of  nonequilibrium.  For  example, 
in  our  experiments  with  plasma  boundary  layers,  significant  electron- 
density  nonequilibrium  can  not  be  adequately  explained  by  use  of  recom¬ 
bination  coefficients  that  assume  only  electrons  as  the  third  body 
for  recombination  (i.e.,  e  +  M'*’  +  e  +  M  +  e,  where  M'*’  and  M  denote  an 
alkali  metal  ion  and  atom,  resp.).  Existing  values  of  heavy-body  recom¬ 
bination  coefficient  Uj  (for  e  +  M"*"  +  Xj  -►  M  +  Xj,  where  Xj  is  an  atom 
or  molecule)  improve  the  agreement  but  remain  unsatisfactory. 

Experimental  studies  of  the  recombination  process  in  which  the 
third  body  is  an  atom  or  molecule  have  been  carried  out  using  (hydro¬ 
carbon)  flames  [1.2],  where  the  products  of  combustion  contained  several 
molecular  species.  To  extract  information  about  the  individual  molecules 


Xj,  the  assumption  was  made  that  the  contributions  from  the  different 
molecules  in  the  mixture  were  additive.  Expressed  equivalently,  in  terms 
of  recombination  coefficients,  it  was  assumed  that  the  overall  rate  a 
(which  was  measured,  in  effect)  could  be  written  in  terms  of  the  indi¬ 
vidual  rates  Oj  by  the  formula 


where  Xj  «  denotes  Che  mole  fraction  of  species  Xj.  Here  nj  is  the 

number  density  of  species  Xj  and  n  is  the  total  number  density. 

The  results  obtained  in  this  fashion  have  been  strongly  criticized 
by  Bates  [1.3],  who  claimed  that  the  recombination  coefficient  for  a 
molecular  mixture  is  not  additive,  but  that  it  is  given  by  the  formula 

®  “  q  ^  ^  Z  X  B  ^p)^  * 

P  j  J  J 

The  quantities  Bj(p)  are  related  to  the  rates  of  de-excitation  of  M  by 
Xj  from  all  energy  levels  above  and  including  p,  to  all  levels  below  p. 
The  factor  r|  is  given  in  terms  of  atomic  constants  and  Che  temperature  T 
as  (2iimekT/h2)3/2  .  kt3/2. 

For  a  single  species,  one  obtains  from  this  expresison  for  a  Che 
result 


a 


j 


It  is  therefore  apparent  Chat  Che  value  of  a  for  a  mixture  of  molecules 
cannot  be  written,  in  general,  in  terms  of  the  individual  Oj.  To  be  able 
to  calculate  a  for  a  mixture,  Che  factors  Bj(p)  must  be  used.  However 
before  being  able  to  rely  on  the  calculated  factors  Bj(p),  it  is  neces¬ 
sary  Co  test  Che  theory  by  obtaining  data  in  which  only  a  single  molecu¬ 
lar  species  Xj  is  present  in  the  recombining  gas  together  with  the 
alkali  metal  M.  It  is  clear  that  flame  experiments  are  unable  to  yield 
such  data. 


Scientific  Merit 


Recombination  and  Ionization  are  fundamental  processes  that  play  a 
role  In  nearly  all  applications  and  natural  phenomena  that  Involve 
partially  Ionized  plasmas.  For  nonequilibrium  situations,  knowledge  of 
the  rates  for  these  processes  are  needed  for  estimating  electron  den¬ 
sity,  and  thereby,  electrical  conductivity.  This  Information  Is  neces¬ 
sary  In  designing  and  understanding  the  performance  of  new  devices,  as 
well  as  In  understanding  the  nature  of  many  terrestrial  and  solar  phe¬ 
nomena. 


Status  Report 

To  achieve  conditions  necessary  for  acquiring  basic  data  and  for 
testing  Bate's  theory,  experiments  are  planned  where  the  basic  approach 
will  be  to  employ  a  supersonic  expansion  of  a  high  temperature  gas. 
The  experimental  apparatus  for  this  work  is  Indicated  schematically 
In  Figure  1.1.  An  arc  heater  will  be  used  to  produce  a  flow  of  high 


Figure  1.1.  Proposed  experimental  aparatus  for  measurement  of 
the  molecular  3-body  recombination  coefficient. 


temperature  Ar  Into  a  large  plenum  chamber.  Potassium  vapor  will  be 
added  to  the  gas,  together  with  various  molecular  species  Introduced 
either  singly,  or  In  controlled  combinations.  The  purpose  of  the  plenum 
Is  to  permit  a  thorough  mixing  of  the  gases  and  to  enable  the  mixture  to 
come  to  thermodynamic  equilibrium. 

The  gas  mixture  then  flows  through  a  chamber  viewport  where  the 
stagnation  temperature  of  the  gas  can  be  measured  with  thermocouples  and 
by  spectroscopic  methods.  A  convergent-divergent  supersonic  nozzle 
follows.  In  which  the  gas  will  be  rapidly  cooled,  so  that  electron 
recombination  Is  Induced.  The  electron  number  density  will  be  measured 
spectroscopically  at  the  exit  of  the  nozzle,  thereby  providing  data  from 
which  the  recombination  coefficient  can  be  calculated.  The  gas  will  then 
pass  through  a  supersonic  diffuser  and  be  exhausted  to  the  atmosphere. 

To  design  the  apparatus.  It  Is  necessary  to  determine  approximately 
the  optimal  experimental  conditions,  within  the  constraints  Imposed  by 
the  operating  limits  of  the  arc  jet.  These  conditions  Include  the  plenum 
pressure  and  temperature,  the  number  densities  of  the  potassium  and  the 
molecular  gas  additives,  and  the  geometry  of  the  converging-diverging 
nozzle.  Since  3-body  recombination  will  occur  In  which  both  the  elec¬ 
trons  and  the  molecules  serve  as  third  bodies.  It  Is  desirable  that 
experiments  be  conducted  for  conditions  such  that  the  latter  process  Is 
dominant. 


At  higher  pressures,  colllslonal  atomic  processes  dominate  radia¬ 
tive  processes,  so  that  one  may  write  Ug 


ngOgCT)  and  aj  =  njajCT), 


where  ag(T)  and  ajCf)  are  functions  of  temperature.  (Here  Og  and  nj 
denote  the  electron  and  molecular  number  densities.)  Shown  in  Figure 
1.2,  are  the  values  of  CjJ  for  N2,  as  calculated  by  Bates,  Malavlya,  and 
Young  [1.4],  and  the  values  of  as  given  by  Hinnov  and  Hirschberg 
[1.5].  We  shall  take  N2  as  the  molecular  third  body  in  the  discussion 
that  follows. 


The  quasi  one-dimensional  steady  state  electron  continuity  equation 
that  governs  the  behavior  of  Og,  is 

dn  .  , 

u  -3 —  +  n  -j—  +  un  T—  An  A  =  n  , 
dx  e  dx  e  dx  e 


-  6  - 


where  denotes  the  net  electron  production  rate,  and  Is  given  by  the 
expression 


n  ■  fn^-S  -  n^l  (n  a'  +  n„  a '  ] 
e  '  iC  e-  '■  e  e  N_  N--* 


2  2 


Here  s*s(T)  *  exp(-  tj^/k-T),  where  is  the  ionization  energy  of 

potassium.  In  obtaining  n^,  the  ionization  rates  have  been  expressed  as 


s(T)  a'  and  S 


N' 


s(T)  by  use  of  the  method  of  detailed 


balancing.  The  gas  flow  speed  u  (as  well  as  all  other  gas  properties)  is 
given  in  terms  of  the  cross-sectional  area  of  the  nozzle  A  =  A(x),  by 
the  well-known  relations  of  gas  dynamics  [1.6]. 


Before  discussing  some  of  the  results  of  integrating  the  continuity 
equation  for  the  ng(x),  it  is  useful  to  introduce  two  concepts  that 
serve  as  a  guide  in  selecting  the  initial  condition  at  the  entrance  to 
the  nozzle.  We  denote  the  ratio  of  the  recombination  rate  by  N2  to  that 


by  electrons  as  r  =  nj^^  ratio  of  the  characteristic 

time  for  recombination  by  N2  to  the  residence  time  of  a  fluid  element  in 


the  nozzle  as  x  =  /‘'res*  ''res®  where 

2  2 1  /  ^  2  2 

L  is  the  nozzle  length  and  a  =  (yRT)*-'*^  is  the  speed  of  sound. 


We  would  anticipate  that  desirable  upstream  gas  states  would 
satisfy  the  conditions  r  ^  1  and  x  ~  1.  (For  x  «  1  the  flow  would 
approach  the  "equilibrium"  limit,  in  which  ionization  and  recombination 
rates  are  approximately  equal.  For  x  »  1  the  flow  would  be  approxi¬ 
mately  "frozen",  where  little  recombination  by  N2  occurs.)  The  condition 
on  r  would  ensure  that  recombination  rate  by  N2  would  exceed  that  by 
electrons. 


From  the  equation  defining  r  we  may  write 


n  = 
e 


and  substituting  for  n  in  the  equation  defining  x,  we  obtain 


[(t)  ,r--  200  /  T  ] 


1/2 


corre- 


(We  have  used  here  the  approximate  expression  -  0.005/  /f~, 

spondlng  to  a  nozzle  length  of  0.1m.)  Shown  In  Figure  1.2  are  the  values 
"N2  ®  function  of  temperature  for  values  of  the  parameter 

(r/x)  equal  to  1  and  10,  and  the  corresponding  values  of  Og  for  r  *  1 
and  10.  Also  shown  are  constant  partial  pressure  lines  for  N2  and  K  with 
Pjj^  =  6  atm  and  Pjj  “  0.1  atm.  These  values  are  estimates  of  the  approxi¬ 
mate  maximum  values  of  the  nitrogen  and  potassium  partial  pressures  that 
can  be  used  with  the  available  experimental  facilities. 
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Figure  1.2.  Temperature  depend¬ 
ence  of  the  3-body  recombination 
coefficients  with  either  N2  or 
electrons  serving  as  the  third 
body 


Figure  1.3.  Domain  of  experi¬ 
mental  conditions  required  for 
determing  the  3-body  recombina¬ 
tion  coefficient  with  N2. 
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The  estimated  range  of  possible  experimental  values  for  n^^  and  the 
temperature  fall  In  the  region  defined  by  T  ^  2000"K  and  by  the  curves 
labeled  (r/x)  -  1  and  p^^^  >  6  atm.  It  appears  that  values  of  (r/t)  up  to 
about  10  can  be  achieved,  which  should  be  adequate. 

To  obtain  a  more  accurate  picture  than  Is  provided  by  the  preceding 
analysis,  a  program  has  been  written  to  Integrate  the  electron  conti¬ 
nuity  equation  numerically.  It  Is  Intended  that  this  program  be  employed 
to  estimate  the  optimal  experimental  conditions  that  should  be  used  for 
determining  a]J^(T).  The  variation  of  n^  with  streamwlse  location  In  the 
nozzle  Is  shown  In  Figures  1.4  and  1.5  for  several  gas  conditions.  The 
nozzle.  In  all  cases,  was  assumed  to  be  conical  with  a  half -angle  of 
3.2"  for  the  converging  part,  and  1.3"  for  the  diverging  part.  The  four 
curves  In  each  figure  show  the  "frozen"  flow,  the  "equilibrium"  flow, 
the  recombining  flow,  and  the  contribution  to  the  recombination  by  just 
the  electron  3-body  process.  Figures  1.4  and  1.5  are  for  stagnation 
temperatures  of  2000"K  and  3000"K  respectively.  With  reference  to  Figure 
1.4a,  Figure  1.4b  shows  the  effect  of  Increasing  n^  ,  and  Figure  1.4c 
shows  the  effect  of  Increasing  n^.  With  reference  to  Figure  1.5a,  Figure 
1.5b  shows  the  effect  of  Increasing  n^j^,  and  Figure  1.5c  shows  the 
effect  of  decreasing  n^. 

The  term  "N2  effect"  for  each  figure  Is  the  fractional  decrease  In 
n^  attributable  to  3-body  recombination  with  N2>  These  results  show  that 
Increasing  the  values  of  n^j^  and  ng  are  generally  favorable  and  Increas¬ 
ing  T  Is  unfavorable  with  respect  to  Increasing  the  fractional  effect  of 
N2.  Work  Is  In  progress  to  explore  systematically  the  variation  of 
electron  density  for  various  gas  conditions  and  nozzle  geometries. 

Future  work  on  the  mathematical  modeling  of  the  recombination 
process  In  a  rapid  expansion  will  include  the  following:  The  model  will 
be  used  to  estimate  the  accuracy  with  which  the  planned  experiments  can 
determine  the  3-body  recombination  coefficient,  and  possible  diagnostic 
methods  for  measuring  Og  will  be  examined.  In  addition,  the  model  will 
be  Improved  to  account  for  skin  friction,  heat  transfer,  and  the  possi¬ 
bility  that  the  electron  temperature  may  differ  from  the  gas  tempera¬ 
ture.  Experimental  work  will  be  undertaken  to  map  out  the  operating 


1.5  Variation  of  electron  density  through  the  nozzle  for  a  stagnation. temperature  of  3000®K 


characteristics  of  the  arc  Jet  for  the  doaaln  of  gas  conditions  indi¬ 
cated  by  the  theory. 

2.2  Placharge  Kffacta:  Plaat-glcctrodc  Interaction 

Introduction 

Insofar  as  practical  applications  of  plasmas  invariably  Involve  the 
transfer  of  current  betmen  electrodes  and  the  plasma,  and  since  the 
performance  of  plasma  devices  is  comnionly  limited  by  mechanisms 
associated  with  this  current  transfer,  a  better  understanding  of  the 
basic  physical  processes  involved  in  the  plasma-electrode  interaction  is 
highly  desirable.  Among  these  limiting  mechanisms  particular  reference 
should  be  made  to  the  excess  voltage  drop  and  power  dissipation  in  the 
plasma-sheath  boundary  layer,  and  the  electrode  erosion  associated  with 
current  concentrations  at  both  cathodes  and  anodes. 

The  work  described  under  this  heading  is  directed  towards  providing 
a  more  complete  understanding  of  the  plasma-electrode  interaction  for 
the  case  of  flowing,  seeded,  thermal  plasmas  at  pressures  of  the  order 
of  one  atmosphere,  and  involves  both  theoretical  and  experimental 
components. 

Consider  a  channel  flow  of  thermal  plasma  which  might  be  either  a 
noble  gas  or  combustion  products  gas  seeded  to  a  suitable  level  (~  1%) 
with  cesium  or  potassium.  We  may  define  the  electric  boundary  layer  as 
that  region  near  the  channel  walls,  composed  of  high  temperature 
insulators  and  electrodes,  where  charged  particle  nonequilibrium  effects 
become  significant  and  even  dominant.  Normally,  the  walls  have  to  be 
operated  at  temperatures  less  than  that  of  the  core  flow  plasma  because 
of  material  limitations. 

On  physical  and  dimensional  grounds  the  region  between  the  core 
plasma  and  the  electrode  surface  can  be  divided  into  a  number  of 
overlapping  layers  of  decreasing  scale.  Proceeding  from  the  core  plasma 
to  the  electrode  these  layers  may  be  listed  as  follows  (see  Figure  2.1). 


(1)  The  gasdynamlc  boundary  layer  of  scale  6  ~  10“ ^m 


(li)  The  ionization  nonequilibriua  (or  aabipolar  diffusion) 
layer  of  scale  the  characteristic  recombination 

length. 

(Hi)  The  sheath  of  scale  equal  to  the  Debye  length  \q, 
typically  10”®m. 

(Iv)  The  Knudsen  layer  where  free  molecule  flow  takes  over 
from  a  continuum  (fluid)  description  and  the  appro¬ 
priate  scale  is  the  mean  free  path  ~  if  ~ 


Figure  2.1.  Potential  distribution  across  plasma  diode. 
1  Equilibrium  Plasma 
II  Ionization  Non-Equilibrium  Layer 
III  Sheath.  (Not  to  scale) 


It  should  be  emphasized  that  the  four  layers  distinguished  above 
are  overlapping  and  no  sharp  boundaries  can  be  distinguished.  Moreover, 
the  scale  lengths  are  nominal  ones  based  on  frees cream  rather  than  local 
conditions.  In  general  Che  scale  lengths  are  ordered  according  to  6  » 


yy  ~  ~  '^i* 


As  Che  current  density  J  is  increased  from  zero,  Che  potential  drop 
at  the  cathode  increases  and  more  ions  than  electrons  are  collected.  The 
reverse  happens  at  the  anode.  However  these  conditions  only  hold  for 


o 

very  snail  currents  J  <  the  random  Ion  current  10  mA./cm  ).  For 
higher  currents,  of  technical  Interest,  the  cathode  must  emit  electrons. 

At  the  anode  it  might  be  expected  that  large  currents  up  to  J  ~  J^, 
the  random  electron  current  (~  10  A/cm^),  could  be  collected  as  a  result 
of  the  reducing  sheath  potential  drop.  However  it  is  commonly  observed 
that  constrictions  in  the  form  of  anode  glow  spots  occur  at  lower 
currents  and  the  sheath  potential  drop  may  change  sign  as  indicated  in 
Figure  2.1.  The  reason  for  this  behavior  is  unexplained  hitherto. 

From  the  foregoing  discussion  it  is  clear  that  the  physical 
processes  involved  are  extremely  complex  and  experimental  data  are  vital 
to  guide  the  theoretical  modeling. 

Scientific  Merit 

In  most  plasma  devices  the  near-electrode  region  and  the  processes 
by  which  current  is  transferred  between  the  plasma  and  the  electrodes  is 
very  poorly  understood,  despite  the  fact  that  the  plasma-electrode 
interaction  is  frequently  the  performance-limiting  process  in  the  use  of 
the  device.  This  program  is  planned  to  shed  new  insight  into  this 
subject  and  could  well  result  in  performance  Improvements  with  respect 
to  voltage  and  power  losses  and  electrode  erosion  effects. 


Status  Report 

(a)  Theoretical  Work 

The  approach  adopted  for  theoretical  modeling  is  to  start  off  with 
the  simplest  relevant  model  and  to  Introduce  the  complexities  of  the 
complete  problem  one  at  a  time.  To  this  end,  the  plasma-sheath  problem 
for  a  plane  electrode  In  contact  with  a  uniform,  equilibrium,  weakly- 
ionized  plasma  has  been  formulated  in  some  generality.  By  assuming  that 
the  plasma  is  Isothermal  and  stationary,  Che  complexities  associated 
with  Che  gasdynamlc  thermal  boundary  layer  and  Joule  heating  are,  for 
Che  present,  eliminated  and  Che  essential  plasma  physics  aspects  of  the 
problem  are  clearly  revealed. 


This  fornulatlon  in  terns  of  Che  conservation  equations  for  the 
electron  and  ion  concentrations  and  nomenta  (Including  Inertia), 
together  with  Poisson's  equation,  give  a  uniform  description  of  the 
ionization  nonequlllbrlua  layer  anu  sheath.  In  both  the  colllonless  and 
collision  dominated  regimes,  without  having  to  resort  to  an  arbitrary 
division  and  matching  of  the  separate  regions.  The  governing  equations 
have  been  non-dlmenslonallzed  and  have  been  programmed  for  numerical 
solution. 

In  addition,  an  analytic  solution  has  been  derived  In  the  quasi- 
neutral  continuum  approxlisaClon  which  describes  the  Ionization 
nonequilibrium  region  exactly,  and  allows  a  number  of  useful  explicit 
results  Co  be  obtained.  This  work  Is  the  subject  of  a  paper  [2.1]  which 
has  recently  been  published.  The  principal  results  are  summarized  below. 

Expressions  for  Che  profiles  of  the  various  parameters  of  Interest 
are  derived  In  normalized  form  as  follows. 

For  Che  plasma  concentration 

n  ■  tanh  y  ; 

For  Che  net  generation  rate  of  electron-ion  pairs: 

~  2  2  - 
n  ■  p  n  ^sech  y  tanh  y  ; 

For  Che  electric  field: 

E  »  [-  sech^  y  +  j]  coth  y  ; 

For  Che  potential: 

^  *  In  tanh  y  -  J  In  sinh  y  ; 


(2.1) 


(2.2) 


(2.3) 


(2.4) 
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For  Che  charged  particle  fluxes: 


f  -  [-  - —  sech  y  -  J]  ;  (2.5a) 

-  [-  2  sech^  y  +  j]  .  (2.5b) 

Here,  Che  nomallzed  distance  Is  y  -  (y/^g);  the  nonnallzed 
concentration  and  generation  rate  are  n  >  (n/i\.),  n  >  (n/Oo,);  the 
nomallzed  potential  Is  f  ■  (e^/lcT);  the  normalized  electric  field  Is 
E  a  (E/E^)  where  Ej^  =  (kT/eX^),  Is  a  characteristic  electric  field;  the 
normalized  current  density  Is  J  -  (J/Jg),  where  =  a„  Ej^  Is  a 

characteristic  current  density;  and  the  normalized  fluxes  are  a 

(rg/n^nPe^^  and  Fj^  a  (F^/n.pjEjj) ,  respectively. 

The  profiles  of  n  and  n,  which  are  Independent  of  current  are  shown 
In  Figure  2.2,  while  Che  potential  profiles  as  a  function  of  current 
density  are  shown  In  Figure  2.3. 

Certain  Important  deductions  can  be  made  concerning  Che  behavior  of 
Che  solution  as  a  function  of  normalized  current  density.  First,  for  a 
cathode  (J  <  0),  It  Is  clear  from  Eq.  2.5a  that  when  -  J  exceeds  the 
relatively  small  value  (pj^/p^),  Che  cathode  must  emit  electrons.  Second, 
for  an  anode  (J  >  0),  It  follows  from  Eqs.  2.3  and  2.4  chat  the 
potential  fall  changes  sign  from  negative  to  positive  for  J  >  1.  Third, 
for  an  anode.  It  follows  from  Eq.  2.5b  Chat  for  J  >  2  Che  anode  must 
emit  Ions. 

In  addition,  this  solution  In  Che  quasl-neuCral,  continuum 
approximation  has  been  used  to  evaluate  Che  magnitude  of  the  neglected 
terms  In  Che  full  equations,  so  that  Che  points  where  space  charge  and 
Ion  and  electron  Inertia  become  significant  can  be  determined.  This 
enables  one  to  derive  conditions  distinguishing  between  Che  cases  of 
collisionless  or  collision  dominated  sheaths. 

The  chief  utility  of  these  results  is  Co  delineate  Che  main  regimes 
of  Che  plasma-electrode  Interaction  and  Co  guide  Che  development  and 
Interpretation  of  results  obtained  by  numerical  solution  of  Che  complete 
equations. 


Figure  2.2.  Profiles  of  nornal- 
Ized  charged  particle  concentra*- 
tlons  and  net  generation  rate. 


Figure  2.3.  Profiles  of  the 
normalized  potential  as  a  func¬ 
tion  of  normalized  current  dens¬ 
ity:  3  <  0,  cathode;  3  >  0, 
anode.  Broken  lines  indicate  the 
asymptotes  corresponding  to  the 
uniform  electric  field  in  the 
distant  plasma 


Experimental  Work 

The  experimental  work  in  progress  is  based  on  the  use  of  a 
versatile  bench-scale  diffusion  burner  capable  of  producing  a  combustion 
plasma  to  3000K  and  seeded  to  1%  by  mass  of  potassium  vapor  to  produce  a 
thermal  plasma  with  a  conductivity  of  up  to  10  mho/m.  This  facility  was 
built  in  connection  with  an  earlier  NSF  grant  for  the  testing  of  plasma 
diagnostic  techniques. 

Currently  a  chimney,  or  flow  channel,  of  internal  cross-section 
1-1/2  in.  X  3-1/2  in.  is  being  designed  to  fit  above  the  flat-flame 
burner  body.  This  has  ceramic  side-walls  to  reduce  heat  losses  and  so 
minimize  the  decrease  in  plasma  temperature  and  conductivity  with  height 
above  the  burner  surface. 
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A  pair  of  cooled  electrodes,  1  cm  diameter,  project  from  the 
ceramic  sidewalls  Into  the  plasma  flow,  to  give  an  electrode  spacing 
adjustable  between  0.5  and  1.5  cm.  Electrodes  of  both  copper  and 
stainless  steel  will  be  used,  and  their  cooling  will  be  adjustable  to 
vary  the  surface  temperature  up  to  within  a  reasonable  safety  margin  of 
the  respective  melting  points. 

When  the  thermal  design  of  the  flow  channel  and  electrodes  has  been 
tested  and  finalized.  It  Is  planned  to  make  electrical  measurements. 
These  will  be  made  by  holding  the  plasma  and  electrode  thermal 
conditions  In  a  steady  state  and  applying  short  (1-10  msec)  voltage 
pulses  between  the  electrodes  at  a  low  repetition  rate.  A  power  supply 
capable  of  delivering  pulses  up  to  10  A  at  400  V  with  variable  pulse 
length  and  repetition  rate  Is  available  for  this  work.  The  purpose  of 
using  a  low  duty  ratio  Is  to  minimize  electrode  damage  which  Is  likely 
to  be  a  problem  under  continuous  energization  at  the  power  levels  of 
interest. 

Initially,  the  primary  diagnostic  technique  will  be  to  observe  the 
voltage-current  characteristic  as  a  function  of  pulse  duration  to 
determine  the  transitions  from  uniform  to  constricted  modes  of  current 
transfer.  Additional  diagnostic  techniques  such  as  voltage  probes  and 
fast  cine  recording  are  also  being  designed  into  the  experimental  set 


2.3  Interaction  of  Discharnes  and  Fluid 


Introduction 


In  several  space  power  and  propulsion  systems  of  potential  Interest 
to  the  Air  Force,  such  as  HPD  thrusters  and  other  MHD  devices,  fluid 
motions  Induced  by  discharges  are  of  critical  Importance  to  both  the 
understanding  and  performance  of  the  system.  In  any  MHD  device  the 
electromagnetic  field  may  affect  the  plasma  flow  field  In  several  ways. 
Among  these  effects  are  the  following: 


(1)  the  effect  of  Joule  heating  on  plasma  temperatures; 


(2)  the  Hartmann  effect; 

(3)  turbulence  suppression;  and 

(4)  the  generation  of  secondary  flows. 

The  first  three  of  these’  effects  have  been  the  subjects  of  previous 
studies  performed  at  HTGL  ([3. 1-3. 4];  these  are  summarized  In  3.5).  The 
fourth  effect,  secondary  flows.  Is  the  area  In  which  our  current  efforts 
are  primarily  focused.  We  are  also  continuing  our  Investigation  In  the 
area  of  turbulence  suppression  by  magnetic  fields. 

The  physical  mechanism  for  secondary  flow  In  a  confined  plasma  flow 
with  a  transverse  magnetic  field  Is  Illustrated  In  Figure  3.1.  For  sim¬ 
plicity  we  Illustrate  a  distribution  of  the  Hall  current,  J^(z),  which 
one  would  expect  from  a  first-order  analysis.  The  magnetic  field  Is 
assumed  constant.  The  resulting  distribution  of  the  J  z  B  Lorentz  force, 
acting  as  a  body  force  against  the  fluid  In  the  channel  cross-plane, 
then  causes  a  secondary  flow  field  which  Is  characterized  by  large-scale 
vortices.  It  Is  believed  that  the  presence  of  this  sort  of  transverse 
flow-field  may  have  profound  Implications  for  the  performance  of  MHD 
devices  [3. 6-3. 7].  The  existence  of  secondary  flow  has  been  Inferred 
from  data  obtained  at  HTGL  on  temperature  profiles,  voltage  drops  and 
heat  transfer  [3.8]. 

The  goals  of  the  present  work  are  to  obtain  direct  measurements  of 
the  secondary  flow  field  under  conditions  in  which  the  current  discharge 
In  the  Hall  direction  Is  controlled,  and  to  correlate  these  results  with 
simultaneous  measurements  of  the  axial  development  of  the  transverse 
electric  field,  which  we  expect  to  be  significantly  affected. 

Scientific  Merit 

Secondary  flow  and  the  Hartmann  effect  are  closely  related,  in  that 
they  are  both  manifestations  of  the  appearance  of  the  Lorentz  force  in 
the  fluid  momentum  equation  for  a  collision-dominated  plasma. 


Figure  3.1.  Physical  mechanism  for  secondary  flow  In  an  MHD  channel. 

For  a  linear  MHD  device  B  Is  assumed  Co  point  In  the  z-dlrectlon;  thus 
Che  LorenCz  force  has  components  In  Che  x-  and  y-dlrectlons. 

The  axial  component,  JyB,  gives  rise  to  the  Hartmann  effect,  while 
the  transverse  component,  “Jx®»  causes  secondary  flow.  Vfhereas  the 
Hartmann  effect  has  been  of  Interest  since  Hartmann's  first  Investiga¬ 
tion,  In  1937,  It  Is  only  recently  that  secondary  flow  has  received 
attention.  The  fact  that  theoretical  models  predict  that  secondary  flow 
may  be  an  Important  effect  In  real  MHD  devices  underscores  the  need  for 


an  experimental  study  of  this  phenomenon.  Except  for  our  own  recent 
work,  no  measurements  of  HHD  secondary  flow  velocities  have  been 
reported. 

Status  Report 

A  preliminary  series  of  experiments  was  conducted  at  HTGL  In  Novem¬ 
ber,  1983,  In  which  measurements  of  the  MHD-lnduced  secondary  flow  field 
were  obtained.  These  results  will  be  reported  at  the  1984  IEEE  Interna¬ 
tional  Conference  on  Plasma  Science,  to  be  held  Hay  14-16  In  St.  Louis 
[3.9],  and  at  the  22nd  Symposium  on  Engineering  Aspects  of  MHD,  to  be 
held  June  26-28  at  Mississippi  State  University  [3.10]. 

The  experiments  were  conducted  at  the  Stanford  M-2  MHD  facility. 
The  channel  was  connected  In  the  Hall  configuration,  with  an  electrode 
pair  at  the  upstream  end  of  the  channel  connected  to  a  pair  at  the 
downstream  end;  the  net  Hall  current  was  controlled  by  means  of  an 
external  power  supply.  Velocities  were  measured  using  laser  Doppler 
anemometry.  A  schematic  of  the  LDA  system  Is  shown  In  Figure  3.2.  Numer¬ 
ous  measurements  of  the  y-dlrected  velocity  component,  v,  were  made  over 
a  grid  covering  one  quadrant  of  the  channel  cross-plane,  and  under  four 
sets  of  boundary  conditions: 

(A)  with  the  Faraday  field  E  open-clrculted,  and  with  the 
magnetic  field  having  positive  polarity; 

(B)  Ey  open-circuited,  magnetic  field  polarity  reversed; 

(C)  Ey  short-circuited,  B-field  positive;  and 

(D)  Ey  short-circuited,  B-field  negative. 

In  each  of  the  four  cases  the  same  flow  rate,  net  Hall  current,  and 
magnetic  field  magnitude  were  maintained,  corresponding  to  an  effective 
MHD  interaction  parameter  of  =  3.1. 

The  results  are  summarized  schematically  in  Figure  3.3.  In  all  four 
cases  the  secondary  flow  field  is  characterized  by  large-scale  vortical 
structures.  The  peak  measured  velocities  range  from  18  m/s  (Case  A)  to 
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Figure  3.2.  Schematic  of  experimental  apparatus. 

35  m/s  (Case  B);  the  latter  corresponds  to  23%  of  the  mean  axial  veloc¬ 
ity,  which  is  far  greater  than  the  1-2Z  peak  secondary-flow  velocities 
characteristic  of  turbulent  flow  through  rectangular  ducts  in  the  ab¬ 
sence  of  MHO  effects.  Figure  3.4  shows  ietailed  results  for  the  mean  y- 
dlrected  velocity  measurements  for  one  of  the  cases  (Case  B). 

If  we  consider  the  simplest  model  for  MHD-induced  secondary  flow, 
as  shown  in  Figure  3.1,  then  we  expect  cases  of  positive  net  Hall 
current  to  lead  to  the  sort  of  flow  fields  found  in  Cases  B  and  C,  in 
which  the  flow  direction  in  the  core  follows  the  direction  of  the  J  x  B 
Lorentz  force,  which  acts  as  a  body  force  on  the  plasma.  The  generation 
of  vorticity  is  a  consequence  of  the  nonuniformity  in  this  force,  as  can 
be  seen  by  inspecting  the  equation  for  the  axial  vorticity,  Q^,  for  the 
case  where  non-MHD  sources  of  vorticity  are  neglected, 

®'^x 
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Ey=O.C.  B  =  + 


Ey=O.C.  B=- 


E„=S.C.  B  =  + 


Ey=S.C.  B^ 


Figure  3.3.  Summary  of  secondary  flow  fields  measured 
for  the  four  cases. 


Since  Jjj  Is  strongly  positive  In  the  core  and  goes  to  zero  at  the  Insu¬ 
lating  sidewalls,  a  first-order  analysis  assumes  that  is  every¬ 

where  positive  over  the  left  half  of  the  cross-plane,  and  hence  that 
follows  the  sign  of  the  B-fleld.  The  results  for  Cases  B  and  C  are  in 
agreement  with  this  expectation. 

The  results  for  cases  A  and  D  do  not  agree  with  the  simple  model. 
In  Case  A  we  observe  the  expected  clockwise  vortex  in  the  core,  but  in 
addition  there  is  a  counter-rotating  vortex  in  the  region  of  the  side¬ 
wall.  In  Case  D  the  direction  of  vortlcity  over  the  entire  field  is 
opposite  to  the  prediction  of  the  first-order  analysis.  The  likely 
hypothesis  for  such  behavior  lies  In  the  possible  presence  of  axial 


Figure  3«4.  The  measured  secondary  flow  field  for 
conditions  of  Case  B. 


velocity  overshoots  In  the  Insulating  wall  boundary  layer.  This  phenome¬ 
non,  which  Is  a  consequence  of  the  Hartmann  effect,  can  lead  to  the 
presence  of  a  region  In  which  dJ^/dz  Is  negative  even  though  the  sign  of 
Is  positive  In  the  core.  This  comes  about  because  of  the  effect  of 
the  u  X  B  field  on  J^,  as  described  by  the  generalized  Ohm's  law.  That 
such  velocity  overshoots  can  occur  In  the  Stanford  channel  was  demon¬ 
strated  by  Rankin  [3.2]. 

It  Is  clear  from  these  preliminary  results  that  secondary  flow  Is 
Indeed  a  significant  phenomenon  In  MHD  devices.  Presently  an  experiment 
is  being  prepared  in  which  this  line  of  investigation  will  be  further 
pursued.  In  addition  to  demonstrating  reproducibility,  we  expect  this 
next  experimental  series  to  add  to  what  we  learned  from  the  results 
discussed  above.  The  following  features  will  be  Incorporated: 


(1)  The  apparatus  will  be  modified  to  allow  scanning  over  the 
entire  cross-plane.  Instead  of  over  just  one  quadrant.  Thus  we 
will  obtain  a  more  comprehensive  picture  of  the  secondary  flow 
field. 

(2)  We  will  measure  u(z)  profiles  In  order  to  see  whether  our  hy¬ 
pothesis  concerning  velocity  overshoots  may  be  valid. 

(3)  Detailed  measurements  will  be  made  of  plasma  voltage  differen¬ 
tials  In  both  the  axial  and  transverse  directions.  This  should 
demonstrate  the  effect  which  secondary  flows  exert  on  the 
distribution  of  the  channel  electric  field. 

As  mentioned  In  the  Introduction  to  this  section,  we  are  also  con¬ 
tinuing  our  Investigation  of  turbulence  suppression.  Although  the 
damping  of  plasma  turbulence  by  a  magnetic  field  has  been  demonstrated 
both  by  Olln,  who  used  a  specially-designed  Pitot  probe  [3.31,  and  by 
Reis,  who  used  laser  Doppler  anemometry  [3.4],  there  remains  an 
uncertainty  as  to  whether  this  turbulence  suppression  causes  a 
laminarlzatlon  of  the  mean  axial  velocity  profile;  i.e.,  u(y)  in  the 
absence  of  current.  Olin's  data  suggest  that  such  an  effect  does  occur; 
Reis's,  that  it  does  not.  As  Che  conditions  for  these  two  investigations 
differed,  we  Intend  to  try  to  resolve  this  question  by  performing  LDA 
measurements  under  conditions  similar  to  those  which  obtained  for  Olln. 
This  should  allow  us  to  determine  whether  his  results  were  due  to 
Inadequacies  In  his  measurement  technique,  or  Co  differences  In  the  flow 
conditions  between  the  two  experiments. 
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